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ABSTRACT 

We study the impact of the spectral dependence of the hnear polarization rotation induced by an 
achromatic half-wave plate on measurements of cosmic microwave background polarization in the 
presence of astrophysical foregrounds. Wc focus on the systematic effects induced on the measure- 
ment of inflationary gravitational waves by uncertainties in the polarization and spectral index of 
Galactic dust. Wc find that for the experimental configuration and noise levels of the balloon-borne 
EBEX experiment, which has three frequency bands centered at 150, 250, and 410 GHz, a crude 
dust subtraction process mitigates systematic effects to below detectable levels for 10% polarized dust 
and tensor to scalar ratio of as low as r = 0.01. We also study the impact of uncertainties in the 
spectral response of the instrument. With a top- hat model of the spectral response for each band, 
characterized by band-center and band-width, and with the same crude dust subtraction process, we 
find that these parameters need to be determined to within 1 and 0.8 GHz at 150 GHz; 9 and 2.0 GHz 
at 250 GHz; and 20 and 14 GHz at 410 GHz, respectively. The approach presented in this paper is 
applicable to other optical elements that exhibit polarization rotation as a function of frequency. 

Subject headings: cosmic microwave background — instrumentation: polarimeters — methods: data 
analysis 



1. INTRODUCTION 

The Cosmic Microwave Background (CMB) polariza- 
tion field can be decomposed into two orthogonal E and 
B modes. On large angular scales the B-mode signal en- 
codes information about infl ation, a period of rapid ex- 
pansion in the early univer se (jKamionkowski et al.lll997l : 
ISeliak fc Zaldarriagalll99"7l) . The signal is characterized 
by the tensor-to-scalar ratio r which quantifies the rela- 
tive strength of inflationary gravitational waves (IGW) 
and density perturbations generated by inflation. The 
level of the IGW signal encodes information about the 
energy scale at which i nflation occurred. Th e current 
upper limit is r < 0.2 (jKomatsu et al.ir201lD . A num- 
ber of experimental efforts are ongoing to search for the 
signal at levels as low as r ~ 0.01 over the coming years. 
On small angular scales, the B-mode signal is dominated 
by the 'lensing signal', which results from gravitational 
lensing of the CMB photons by the large scale structure 
of the univer se. The lensing converts E-m odc to B-mode 
polarization (iZaldarriaga fc Seliaklll998l ). 

Galactic foregrounds are expected to be a source of 
confusion for measurements of the B-mode signal. Above 
70 GHz the polarized emission from Galactic dust is pre- 
dicted to dominate over much of the sky and be compa- 
rable to the IGW signal with an r va lue of 0.1 or less 
even for th e cleane st regions of the sk y ([Page et al.ll2007l : 
iGold et al.i 2009: Fraisse et al.ll2011[ ). Therefore many 
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experimental efforts plan to employ multiple frequencies 
which will enable foreground identification and subtrac- 
tion. 

Some CMB polarimeters use a half- wave plate (HWP) 
to modula te the observed linear polariza tion, such 
as EBEX ( Rcichb orn-Kiennerud e t al.l l2010D. SPIDER 
(iRunvan et al.. .20101 ) and POLARBEAR ()Arnold et al.l 
|2010( ). When observing at multiple frequency bands si- 
multaneously, an achromatic half- wave plate (AHWP) 
can be used. An AHWP is a stack of monochro- 
matic HWPs with a part icular set of orienta, tion an- 
gles relative to each other ()Pancharatnam|[l955[ ) . While 
an AHWP has a higher modulation efficiency across a 
broad frequency range compared to a single HWP, it 
rotates the polarization angle of the inc ident light by 
an a, r nount that de pends on frequency (jHananv et al.l 
120051 : iMatsumura et al. 200^. The amount of rotation 
depends on the construction parameters of the AHWP 
and on the spectrum and polarized intensity of the con- 
stituent signals, which for this paper are CMB and Galac- 
tic dust. With knowledge of the spectrum and relative 
polarization intensity, the amount of rotation can be cal- 
culated and corrected. However, while the spectrum of 
the CMB component is well known, that of dust is not. 
The polarized intensities of dust and CMB are also not 
well known. These uncertainties may pose challenges in 
the extraction of the underlying IGW signal. Various 
authors studied the impact of HW P non-ide alities on 
measurements of CMB polarization (jBrown ei al. 2001 
iBrvan et al.ll2010f ). However, this particular frequency 
dependent rotation effect has not been studied in the 
context of B-mode measurements. The goal of this paper 
is to quantitatively assess this effect. For concreteness 
we adopt the AHWP model, frequency bands and ap- 
proximate noise information that are applicable to the E 
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and B experiment (EBEX) ()Reichborn-Kiennerud et "aLl 
l201Clf ). a balloon-borne CMB polarimeter targeting the 
IGW signal at the r - 0.04 level. 

In Section 2 we describe the basic components of the 
simulation. Section 3 focuses on quantifying the effect 
of rotation due to the AHWP in the 150 GHz band. In 
Section 4 we use multiple frequency information to ac- 
count for rotation due to Galactic dust. In Section 5, we 
study the additional effects of uncertainties in the spec- 
tral response of the instrument, and in Section 6 we make 
concluding remarks. 

2. DESCRIPTION OF THE SIMULATION 

We simulate input Stokes Q and U signals due to the 
CMB and Galactic dust emission on a 10° x 10° area of 
the sky centered on {I, b) = (252°, —52°) in Galactic coor- 
dinates which is close to the center of the area targeted by 
EBEX. The maps are smoothed with an 8' FWHM Gaus- 
sian beam then projected to a fiat sky and pixelized with 
a square 6.9' pixel. Same simulations with a 20° x 20° 
patch in the same region validate that conclusions pre- 
sented in this paper do not depend on patch size. The 
input CMB polarization angular power spectra, includ- 
ing both the prim ordial and lensing signal, are gener- 
ated with CAMB (jLewis et al.ll2000li us ing the best fit 
WMA P 7- year cosmological parameters (jKomatsu et al.l 
120 lit) and r = 0.05, unless otherwise indicated. Our po- 
larized f oreground simulatio n follows the prescription de- 
tailed in lStivoli et al.l (|2010f ) and is briefly reviewed here. 
The dust in tensity and its frequency scaling are given by 
'model 8' of lFinkbeiner et ahl (|1999[ ). The dust polariza- 
tion fraction is modeled for cases of 2%, 5%, and 10%. 
A polarization fraction higher than 10% would exceed 
the limit based on WMAP observati ons at interm ediate 
and h igh Galactic latitudes (jPage et a l. 2007; Gol d et al.l 
|2009| ). Both the dust polarization fraction and the fre- 
quency scaling are assumed to be uniform over the sim- 
ulated sky area. Ob servations suggest that this is a 
good approxim ation f Planck Collaboration et ahl 120111 : 
iFinkbeiner et a l. 1999). The pattern of the polarization 
angles on large angular scales (/ < 1 00) is given by th e 
WMAP dust polarization template (jPage et al.l l2007f ). 
On smaller angular scales (/ > 100) we add a Gaussian 
fluctuation power ad opting a recipe first presented by 
iGiardino et all (|2002D . Figure [1] shows the power spec- 
tra of the CMB and of Galactic dust. For a level of 5% 
fractional polarization the expected level of Galactic dust 
is comparable to the B-mode signal at £ = 90. 

To simulate the operation of the AHWP we 
use the Mueller matrix formalism as described by 
iMatsumura et all ()2009D . The level of input polarized 
signal is calculated for each map pixel in 50 frequency 
bins for each of the experiment's three top- hat bands 
(see Table [Ij . For each map pixel the detected intensity 
as a function of AHWP angle, which we call intensity vs. 
angle (IVA), is calculated for each frequency bin with an 
angular resolution of 0.05° and the total per-band IVA 
is the average of the 50 IVAs. The detected polarization 
angle, which is rotated relative to the input polarization 
angle, is encoded by the phase of the band-averaged IVA. 
To obtain the rotated map observed by the detector, we 
multiply each pixel of the input maps by a rotation ma- 
trix with the calculated phase of the band-averaged IVA. 
The frequency and IVA angular resolution are chosen to 
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Fig. 1. — CMB (solid) and Galactic dust (dashed) B-mode power 
spectra at 150 GHz assuming an r value of 0.05. The spectra from 
dust are for the specific area of sky simulated in this work and are 
given for three fractional polarization cases of 2, 5, and 10 %. 



optimize computation time while giving negligible bias in 
the results. The construction parameters of the AHWP 
are given in Table [T] 



Indices of refraction of AHWP 


rio = 3.047, rie = 3.364 


Thickness of each wave plate 


1.69 mm 


150 GHz band 


133 - 173 GHz 


250 GHz band 


217 - 288 GHz 


410 GHz band 


366 - 450 GHz 


Band shape 


top- hat 


Orientation angles of 5-stack AHWP 


(0°, 25°, 88.5°, 25°, 0°) 



TABLE 1 

AHWP AND BAND PARAMETERS USED IN THE SIMULATIONS. 



We calculate both EE and BB power spectra simulta- 
neously using the flat-sky approximation (JKaisor 1992). 
Each simulation is run 100 times with different CMB and 
noise realization, unless otherwise noted. In this study 
we focus on the BB power spectra. The result quoted 
for a given £ bin is the mean of the 100 simulations and 
the error bar is the standard deviation. Figure [2] shows 
a validation of the process of generating CMB Q and 
U maps and estimating the underlying E and B-mode 
power spectra. No rotation due to the AHWP has been 
included in this validation. 

For simulations that include the effects of instrumen- 
tal noise we assume it is homogeneous and has a white 
spectrum, and add its realization to the signal to make 
a combined input map. In our simulation we use an in- 
strumental noise per pixel of 1, 2.8, and 25 /xKcmb for 
the 150, 250, and 410 GHz bands, respectively. Figure [2] 
shows a validation of the noise generation and estimation 
process. 

3. THE EFFECT OF GALACTIC DUST 

If the shape of the frequency band is known then the 
rotation induced on the CMB signal alone can be calcu- 
lated and compensated exactly because the spectrum of 
the CMB is known. This rotation is uniform across the 
sky and with the parameters given in Table [T] is 55° in 
the 150 GHz band. The presence of Galactic dust mod- 
ifies the intensity and angle of the net incident polar- 
ization and thus the amount of rotation induced by the 
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Fig. 2. — Validation of signal and noise power spectrum estima- 
tion. We generate 100 Q and U maps using an underlying power 
spectra (black solid curves) and use a flat-sky approximation to 
estimate the power spectra (red triangles with error bars). The 
size of error bars agrees with predictions for the contributions of 
cosmic and sample variance. We also make 100 noise only Q and U 
realizations using white noise with RMS of 1 /iK (magenta dashed 
line) and estimate the power spectrum (blue solid circles with error 
bars). The estimated noise power spectrum is shown slightly offset 
in I to enhance clarity. 

AHWP. The spectral dependence and spatial distribu- 
tion of Galactic dust polarization is not precisely known 
and therefore the amount of rotation it induces can only 
be estimated. How big is this extra rotation? Can it sim- 
ply be ignored because it is negligible? In the remainder 
of this section we assess these questions for the 150 GHz 
band. 

These first simulations include CMB and dust, without 
instrumental noise. We calculate rotated Q and U maps 
resulting from passing the total (CMB+dust) incoming 
polarization through the AHWP. We then 'de-rotate' the 
maps by the calculated rotation angle for CMB only, sim- 
ulating ignorance of the effects of the dust foreground on 
the rotation. We subtract the input dust Q and U maps 
from the de-rotated map to acknowledge the presence of 
the effects of dust polarized intensity on the total Q and 
U maps. Note however that the effect of rotation due to 
dust polarization, which is a consequence of the AHWP, 
is left in the map. We then calculate the angular power 
spectrum of the resulting maps for 2%, 5%, and 10% of 
dust polarization (see Figure [3]). For the fiducial value of 
r = 0.05 ignoring the effect of rotation introduces notice- 
able bias in the estimation of the CMB power spectrum 
for 10% dust polarization but not for 2% dust polariza- 
tion. For 5% dust polarization the bias is only noticeable 
in the lowest two and the highest £ bins. 

4. REMOVING AHWP INDUCED ROTATION IN 
DUST SUBTRACTION 

In the previous section we showed that for levels of 
polarized dust of more than 5% the effect of the rotation 
due to dust in the AHWP cannot be ignored. In this 
section we employ a simple form of dust subtraction in 
an attempt to correct for the rotation. In this approach 
we make two initial assumptions in order to extract the 
dust frequency scaling information: 

1. the signal at the 150 GHz band is dominated by 
the CMB and dust can be neglected; 

2. the signal at the 410 GHz band comes entirely from 
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Fig. 3. — Effect of rotation due to 2%, 5%, and 10% polarized 
dust on the estimation of the B-mode power spectrum. The in- 
put CMB power spectrum (red solid circles) follows the underlying 
assumed power spectrum (black solid line). The difference be- 
tween the input CMB and the power spectrum of the map after 
de-rotation by the rotation angle corresponding to the CMB alone 
(magenta triangles) indicates the effect of polarized dust. After 
subtracting the dust intensity, only the rotation due to the pres- 
ence of dust remains (blue diamonds). The power spectrum of a 
map of the difference between the input map and the dc-rotatcd, 
dust-subtracted map (yellow crosses) quantifies the effect of the 
rotation due to dust alone. Power spectra points are all calculated 
at the same I bins but are shown slightly offset in I to enhance 
clarity. To reduce clutter we only plot data for every other I bin 
for I > 120. 



dust and CMB can be neglected; 

We prepare the total (CMB-|-dust) rotated polariza- 
tion maps in 150, 250, and 410 GHz bands and add noise 
in the map domain. We then calculate the polarization 
intensity maps and the signal RMS in all three bands 
assuming that the noise RMS is exactly known. Follow- 
ing assumption 1 the map RMS of CMB at 150 GHz is 
known. We extrapolate the CMB level to 250 GHz and 
calculate the map RMS for dust at 250 GHz. Following 
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assumption 2 we also obtain the map RMS for dust at 
410 GHz. Using the dust levels at 250 and 410 GHz we fit 
a gray body dust model, given by a power law multiplied 
by an 18 K blackbody. The power law spectral index 
is taken to be uniform across the entire simulated sky 
area. The top hat spectral response of the instrument 
is assumed to be precisely known (we address uncertain- 
ties in the bandpass in Section [5]). The fitted dust model 
is used to calculate the level of dust at the 150 GHz 
band, extrapolated from the 410 GHz map, and to cal- 
culate and correct for the rotation angle at this band due 
to the combination of dust and CMB. Wc make a final 
map that contains an estimate of the CMB alone after 
corrections for both dust polarized intensity and rotation 
induced by the AHWP. We calculate the power spectrum 
of this map including subtraction of an estimate of the 
noise spectrum. For an estimate of the noise spectrum 
we use the known input RMS. As a test of the entire 
pipeline we run it with no dust and no noise and vali- 
date that the extracted power spectrum agrees with the 
input CMB power spectrum. 

For each set of 100 simulations, we determine whether 
the final estimated CMB power spectrum is biased or not. 
The power spectrum is conservatively assumed biased 
if the mean power estimated in any £ bin is outside of 
the l-cr cosmic variance error bar. We find that for the 
nominal noise levels, r = 0.01 or above, and all dust 
polarized fractions at or below 10% the dust subtraction 
procedure recovers an unbiased estimate of the B-modc 
power spectrum. The results for 10% dust polarization 
fraction are shown in Figure 21 Only data points with 
signal to noise ratio (SNR) > 1 are plotted. For r = 0.009 
and lower, and 10% dust polarization fraction, we find 
that the recovered B-mode power spectrum is biased at 
the lowest £ bin. When the dust polarization is lower 
the r level that can be recovered without bias is higher 
because the higher relative noise at 410 GHz has a larger 
effect on the CMB estimate at 150 GHz. For 2% dust 
polarization fraction, we can recover r as low as 0.02. 

We also use a different approach to quantify the bias 
caused by the dust subtraction procedure. Wc run 100 
simulations with an input of r = and 10% dust po- 
larization fraction. We fit a non-zero r to the difference 
between the estimated and input CMB power spectrum 
at the lowest £ bin, while keeping the shape of the pri- 
mordial B mode signal. We consider this fit as the lower 
limit of r value we can detect using this dust subtraction 
method. We find a best fit with r = 0.01, which is close 
to the result we found earlier. 

The error bars shown in Figure U] include only cosmic 
variance for the input CMB, but include the total error 
on the recovered CMB power spectrum (including the 
effect of correcting for the rotation). On average the 
total error is larger by 30% compared to the error from 
the combination of cosmic variance and instrument noise 
(Fig. 131 upper panel) . 

5. UNCERTAINTY IN DETECTION BAND AND 
HIGH FREQUENCY SPECTRAL RESPONSE 

So far we assumed that the spectral response of the in- 
strument is known. Only the frequency scaling of Galac- 
tic dust is determined from the fit. Uncertainty in the 
spectral response leads to uncertainty in the amount 
of rotation induced by the AHWP. To assess the level 
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Fig. 4. — Comparison between the underlying CMB model (black 
line), the input CMB (red circles, error bars include only cosmic 
variance) and the estimate of the CMB power spectrum using a 
map in which the effects of both dust polarized intensity and dust 
induced extra rotation has been accounted for (blue diamonds). 
The final error bars are on average 30% larger than those due to 
cosmic variance and instrument noise (magenta dot-dashod error 
bars, only plotted in three mid-^ range bins for clarity). With 10% 
dust polarization, a crude dust subtraction algorithm (see text) can 
account for the rotation induced in the AHWP with an r value as 
low as 0.01 (top panel). In the bottom panel the input r is zero and 
the underlying CMB spectrum has only a Icnsing signal; Galactic 
dust is 10% polarized. The estimated CMB spectrum has a best 
fit r = 0.01. To reduce clutter we only plot data at every other £ 
bin for £ > 800. 



of this effect quantitatively we assume a top-hat band 
shape that is characterized by two parameters, center 
and width. We simulate CMB and dust signals with 
nominal bands, and analyze the maps using the dust sub- 
traction algorithm discussed in the previous section but 
assuming varying band-widths, or varying band-centers 
(but not both simultaneously). No instrumental noise is 
included. All simulations have 10% polarized dust and 
r = 0.05. Wc search for the level of shift in band-center 
or change in width that leads to bias in the estimation 
of the final CMB spectrum. We use the same criterion 
for bias as described in Section [H 

5.1. Shift of Band-Center 

Simulations are carried out by shifting only one band- 
center, keeping the other two fixed at their nominal val- 
ues. We find that shifts of more than 1, 9, and 20 GHz 
for the 150, 250, and 410 GHz bands, respectively, lead to 
biased power spectra. The limit for shift of the 150 GHz 
band is due to mixing between E and B modes: an er- 
ror in band-center leaves the CMB slightly rotated after 
correction for the AHWP is applied (using the nominal 
band) and thus a portion of the E-mode signal is mixed 
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into the B-mode signal. This is apparent in Figure [S] 
(panel b), which shows a 2 GHz shift for the 150 GHz; 
the bias is primarily at high t. The limits on the 250 and 
410 GHz bands mainly come from misestimate of rota- 
tion due to dust, but because dust is not dominant at the 
150 GHz band the requirement is less stringent. Panels 
(c) and (d) in Figure [S] show shifts of 10 and 22 GHz for 
the 250 and 410 GHz bands, in which bias due to dust is 
found only at the lowest i bin. 

5.2. Misestimate o] Band-Width 

Simulations are carried out by changing one band- 
width at a time, keeping the other two fixed at their 
nominal values. We find that a misestimate of band- 
width by more than 0.8, 2, or 14 GHz for the 150, 250, 
and 410 GHz band, respectively, exceeds the criterion for 
no bias. The result for the 150 GHz band with a change 
of 1 GHz in width is shown in Figure [6l The cause of 
bias in any of the bands is a misestimate of total power 
detected at the particular band thus a misestimate of 
both the polarized signal intensity and the rotation due 
to dust. For this reason the bias is largest at the lowest 
I bins. 

5.3. Effects of High Frequency Spectral Leak 

Dust intensity is rising up to ~ 2 THz. A higher than 
expected and unknown instrumental response at out-of- 
band frequencies, which is called a 'spectral leak', may 
bias the subtraction of the dust signal and by extension 
the estimate of the underlying CMB signal. We simu- 
late two specific leaks, which are both top-hat in shape, 
a narrow leak between 1750 and 1850 GHz and a broad 
leak between 500 and 2000 GHz. For both cases the 
power in the leak is adjusted to be 0.1%, 1% and 1% 
of the in-band power for the 150, 250, and 410 GHz 
bands, respectively, as measured with a 300 K black- 
body source. These values are readily achievable ex- 
perimentally (jPolsgrovell2009f ). We properly include the 
change i n the refraction indices of sapphire with fre- 
quency (jLoewenstein et all |1973() . Maps are prepared 
with signals that include power in the leak, but are ana- 
lyzed, including the steps of dust subtraction, assuming 
no knowledge of the leak. Instrumental noise is not in- 
cluded in the simulation. For both cases we find no biases 
in the estimate of the final CMB power spectrum. Figure 
[7] shows the case for the broad leak. 

6. DISCUSSION AND SUMMARY 

The spectral response of an achromatic half- wave plate 
may induce bias in the estimation of polarized signals. 
We analyze the level of such bias in the context of mea- 
surements of the B-mode signal of the CMB in the pres- 
ence of Galactic dust, the dominant source of foreground 
emission in cases of interest here. For concreteness, we 
use the specific experimental configuration correspond- 
ing to the EBEX balloon-borne experiment. 

For the area of sky considered we find that with rea- 
sonable assumptions about the magnitude and spectral 
shape of dust, the effects of rotation induced by the 
AHWP are only appreciable when dust is polarized at 
a level of about 5% and above and the tensor-to-scalar 
ratio r is less than ~ 0.05. In the regime when the ef- 
fects of rotation are appreciable, even a crude process of 
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Fig. 5.— Effect of band-center shift of 1 GHz of the 150 GHz 
band (panel a), 2 GHz of the 150 GHz band (panel b), 10 GHz 
of the 250 GHz band (panel c) and 22 GHz of the 410 GHz band 
(panel d), respectively. In panel (b) excess power at high £ comes 
from mixing of E and B-mode. In panels (c) and (d) it comes from 
misestimate of the effects of dust. Figure symbols are the same as 
in Figure in 
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Fig. 6. — A 1 GHz increase in the assumed band- width of the 
150 GHz band relative to the nominal width leads to an overesti- 
mate of the level of dust in this band and thus to a misestimate of 
the rotation due to dust. This leads to excess power at low £. The 
misestimate of the band-width also results in an underestimate of 
the CMB power which leads to a small (less than l-cr) power deficit 
at high i bins. Figure symbols arc the same as in Figure |4] 
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Fig. 7. — The effect of a broad high frequency spectral leak (see 
text for details) is negligible with r = 0.05 and 10% dust polariza- 
tion. Figure symbols arc the same as in Figure U 

dust estimation and subtraction mitigates the effects of 
AHWP rotation to below detectable levels. For exam- 
ple, using the crude dust subtraction process we find no 



bias in the estimation of the B-mode power spectrum for 
dust polarization fraction as large as 10% and r as low 
as 0.01. For 2% dust polarization fraction, r of 0.02 or 
higher is recovered without bias. We also find that the 
dust subtraction causes the power spectrum error bars 
to increase by a modest 30% on average. 

Employing the same dust estimation and subtraction 
process, but now assuming errors in knowledge of the 
experiment's detection band-center and band-width, we 
find the accuracy with which these need to be measured. 
For example, for the particular experimental configura- 
tion considered, we found that band-center and band- 
width of the 150 GHz band need to be determined to 
better than 1 and 0.8 GHz, respectively. It is possible 
that this requirement may not need to be as stringent if 
a more sophisticated foreground estimation and subtrac- 
tion process is used. This research is ongoing. 

We explore the sensitivity of the particular experimen- 
tal configuration to high frequency spectral leaks. Using 
a rejection level that is readily achievable experimentally 
we show that spectral leaks are not expected to pose 
challenges for the operation with an AHWP. 

The analysis and subtraction approach discussed in 
this paper are applicable to other optical elements for 
which polariza tion rotat i on is a function of frequency. 
For example, lO'BrientI ()2010[) describes a broadband, 
mm-wave detection technique that is based on sinuous 
antenna. It is well documented that such antennas 
change the phase response of polarized signals, and that 
this effect is frequency dependent. Thus they exhibit fun- 
damentally the same behavior as an AHWP. Our meth- 
ods and approach apply to such cases. 
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